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Abstract. We define and investigate a notion of entropy for quantum error 
correcting codes. The entropy of a code for a given quantum channel has a 
number of equivalent realisations, such as through the coefficients associated 
with the Knill-Laflamme conditions and the entropy exchange computed with 
respect to any initial state supported on the code. In general the entropy of 
a code can be viewed as a measure of how close it is to the minimal entropy 
case, which is given by unitarily correctable codes (including decoherence-free 
subspaces), or the maximal entropy case, which from dynamical Choi matrix 
considerations corresponds to non-degenerate codes. We consider several ex- 
amples, including a detailed analysis in the case of binary unitary channels, and 
we discuss an extension of the entropy to operator quantum error correcting 
subsystem codes. 



1. Introduction 

Quantum error correcting codes are a central weapon in the battle to overcome 
the effects of environmental noise associated with attempts to control quantum 
mechanical systems as they evolve in time [HE]. It is thus important to develop 
techniques that assist in determining whether one code is better than another for a 
given noise model. In this paper we make a contribution to this study by introducing 
a notion of entropy for quantum error correcting codes. 

No single quantity can be expected to hold all information on a code, its entropy 
included. Nevertheless, the entropy of a code is one way in which the amount of 
effort required to recover a code can be quantified. In the extremal code has 

zero entropy if and only if it can be recovered with a single unitary operation. This 
is the simplest of all correction operations in that a measurement is not required 
as part of the correction process. These codes have been recently coined unitarily 
correctable [SI HI [5], and include decoherence-free subspaces [H |7l|8l|9l [10] in the case 
that recovery is the trivial identity operation. Thus more generally the entropy can 
be regarded as a measure of how close a code is to being unitarily correctable, or 
decoherence-free in some cases. 

In the next section we introduce the entropy of a code, along with required 
nomenclature. We also consider an example motivated by the stabilizer formalism 
[TT] and discuss an extension of code entropy to operator quantum error correcting 
subsystem codes [121 [Mj [15] . We then consider in detail an illustrative class of 
quantum operations for which the code structure has recently been characterised, 
the class of binary unitary channels [TB I [TT l [TS l fT9 l [20 ]. [21 ] [22] . 



2. Entropy of a Quantum Error Correction Code 

Let p denote a quantum state: a Hermitian, positive operator, satisfying the 
trace normalisation condition Trp = 1. A linear quantum operation (or channel) $, 
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which sends a density operator p of size N into its image p' of the same size may 
be described in the Choi-Kraus form [231 El] 

M 

(1) p' = '^{p) = Y^E,pEl 

i=l 

The Kraus operators Ei can be chosen to be orthogonal {Ei\Ej) = TtEJEj = diSij, 
so that the non-negative weights di become eigenvalues of the dynamical (Choi) 
matrix, D$ = {{E.i\Ej)). We refer to the rank of the Choi matrix as the Choi rank 
of <&. Observe that the Choi rank of $ is equal to the minimal number of Kraus 
operators required to describe $ as in 

Hence in the canonical form the number AI of non-zero Kraus operators does 
not exceed iV^. Due to the theorem of Choi the map $ is completely positive (CP) 
if and only if the corresponding dynamical matrix is positive if and only if $ has 
a form as in ([1]). The map $ is trace preserving, Tvp' = Trp = 1, if and only if 

'^f^i EjEi — 1 where we assume some Ei are zero if M is less than N"^. The family 
of operators Ei is not unique. However, if {Fj} is another family of operators that 
determine $ as in ([1]), then there is a scalar unitary matrix U — (uij) such that 
^ij^j f*^^ refer to this as the unitary invariance of Choi-Kraus 

decompositions. 

2.1. Entropy Exchange and Lindblad Theorem. To characterise the informa- 
tion missing in a quantum state one uses its von Neumann entropy, 

(2) S{p) = -Trplogp. 

We will use the convention that log refers to logarithm base two as this provides a 
cleaner operational qubit definition in the context of quantum information. 

In order to describe the action of a CP map <&, represented in the canonical 
Choi-Kraus form ([T]), for an initial state p we may compare its entropy with the 
entropy of the image S{p') ~ S{^{p)). To obtain a bound for such an entropy 
change we can define an operator a = cr(<&, p) acting on an extended Hilbert space 

(3) ay = TvpEjE^, i, j ^ 1, . . . , . 

If the map is stochastic, the operator a is positive definite and normalised so 
it represents a density operator in its own right, a G A^jv^ (specifically, it is an 
initially pure environment evolved by the unitary dilation of $). Observe that for 
any unitary map, ^u{p) — U pU\ the form |T]) consists of a single term only. Hence 
in this case the operator a reduces to a single number equal to unity and its entropy 
vanishes, S'(cr($c/, p)) = 0. 

The auxiliary state a acting in an extended Hilbert space was used by Lindblad 
to derive bounds for the entropy of an image p' = $(/)) of any initial state. The 
bounds of Lindblad [551 , 

(4) < \Sip')-Sia)\ < Sip) < Sia) + Sip'), 

are obtained by defining another density matrix in the composite Hilbert space 

M M 

(5) a; = , 

1=1 j=i 
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where Ad — N'^ and \i) forms an orthonormal basis in T-Lm- Thus, uj is simply 
the system and an initiaUy pure environment evolved by the unitary dilation of $. 
Computing partial traces one finds that Tv^uj = a and Tr^/w = p' ■ It is possible 
to verify that <S'(w) = S'(p), so making use of the subadditivity of entropy and the 
triangle inequality [26] one arrives at (|4|). 

If the initial state is pure, that is if S{p) = 0, we find from (jl]) that the final state 
p' has entropy S{a). For this reason S{'j) was called the entropy exchange of the 
operation by Shumacher :27j. In that work an alternative formula for the entropy 
exchange was proven, 

(6) S{a{'^,p)) = 5((<i>® , 

where jV') is an arbitrary purification of the mixed state, Trsj?/;) (■01 = p. Thus, the 
entropy exchange is invariant under purification of the initial state and remains a 
function only of the initial density operator p and the map $. 

2.2. Quantum Error Correcting Codes. A quantum operation $ allows for an 
error correction scheme in the standard framework for quantum error correction 
[TTl [5S1 [551 [301 [31] if there exists a subspace C such that for some set of complex 
scalars A = (Ay ) the corresponding projection operator Pq satisfies 

(7) PcEJe.Pc = X.,jPc for aU i,j = 1,...,N\ 

Specifically, this is equivalent to the existence of a quantum recovery operation ^ 
such that 

(8) ^o^oVc^Vc, 

where Vc is the map Vcip) — PcpPc- The subspace related to Pc determines a 
quantum error correcting code (QECC) for the map $. A special class of codes are 
the unitarily correctable codes (UCC), which are characterised by the existence of 
a unitary recovery operation '^{p) — UpW. These codes include decoherence-free 
subspaces (DFS) in the case that "if is the identity map, '^{p) — p- 

It can be shown that the matrix A = (Ay ) is Hermitian and positive, and is 
in fact a density matrix, so it can be considered as an auxiliary state acting in an 
extended Hilbert space of size at most N"^ . It is easy to obtain a more refined global 
upper bound on the rank of A in terms of the map $. 

Lemma 1. Let A be the matrix determined by a code C for a quantum map $. 
Then the rank of A is bounded above by the Choi rank of^. 

Proof. Without loss of generality assume the Choi matrix is diagonal. We 
have for all i, 

(9) A,, dimC = Tr(A,,Pc) = Tr{PcEjE,Pc) < TviEjE,) = {E,\E,), 

and the result follows from the positivity of Z?$ and A. ■ 

Unitarily correctable codes are typically highly degenerate codes, as the map $ 
collapses to a single unitary operation when restricted to the code subspace. In par- 
ticular, the unitary invariance of Choi-Kraus representations implies the restricted 
operators EiPc are all scalar multiples of a single unitary. More generally, one can 
see that a code C is degenerate for precisely when the Choi rank of $ is strictly 
greater than the rank of A. Indeed, the Choi rank counts the minimal number of 
Kraus operators required to implement $ via (2.1), and satisfaction of this strict 
inequality means there is redundancy in the description of $ o 7-"^ by the operators 
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EiPc- Thus, for these reasons we shall refer to codes as non- degenerate if the Choi 
rank of <& coincides with the rank of A and if the spectrum of A is equally balanced 
- that is to say that non-degenerate codes correspond to the maximally degenerate 
error correction matrix, A proportional to the identity matrix 

2.3. Entropy of a Code. Assume now that an error correcting code C exists and 
all conditions ([7]) are satisfied. If a quantum state p is supported on the code then 
PcpPc — P and calculation of the entropy exchange ([3]) simplifies, 

(10) (T„ = l^rpElEj = TiPcpPcEjE, = TrpPcEjEjPc = TrpA„Pc = . 

In this way we have shown that the error correction matrix A is equal to the auxiliary 
matrix a of Lindblad, provided the initial state belongs to the code subspace. 

From another direction, given an error correcting code C for a map in [4l I32j 
it was shown that there is a quantum state r and an isometry V such that for all 
P = PcpPc, 

(11) ^{p) ^V{T(g) p)V^. 

The result, which can be seen as a consequence of the decoupling condition of [33], 
gives an explicit way to "see" a code at the output stage of a quantum process 
for which the code is correctable. The result (and its subsystem generalization - 
see below) may also be viewed as a formalisation of the subsystem principle for 
preserving quantum information [7J. From the proof of this result one can see 
directly that the entropy of t satisfies S{t) — S{A). This equality follows also from 
the fact that r and A can be interpreted as the states obtained by partial trace of 
an initially pure state with respect to two different subsystems. 

Thus, from multiple perspectives we find motivation for the following: 

Definition 2. Given a quantum operation <I> with Kraus operators {£"«} and a code 
C with matrix A given by ([7]), we call the von Neumann entropy S'($,C) :— S{A) 
the entropy of C relative to $. 

The entropy of a code depends only on the map and the subspace defined by Pc , 
not on any particular state in the code subspace. Thus, the entropy exchange will 
be the same for all initial states supported on the code subspace and is therefore a 
property of the code itself. 

In the following result we determine what possible values the code entropy can 
take, and we derive a characterisation of the extremal cases in terms of both the 
code and the map. 

Theorem 3. Let ^ be a quantum operation and let C be a code with matrix A given 
by Then S'($,C) belongs to the closed interval [O,log£'], where D is the Choi 
rank of^. Furthermore, the extremal cases are characterised as follows: 
(i) S{^,C) — if and only if C is a unitarily correctable code for $. 
(m) S{^,C) — logD if and only if C is a non-degenerate code for <I>. 

Proof. By Lemma [T] and the subsequent discussion, the maximal entropy case 
occurs when the rank of A and D$ coincide and the spectrum of A is equally 
balanced; that is, the code is non-degenerate. This occurs (by a standard spectral 
majorization argument) precisely when the code entropy satisfies S'($,C) = log 13. 

For the minimal entropy case, first suppose that C is a UCC for Then by ([8|) 
there is a unitary operation U (p) = UpU^ such that $ o —UoVq. Thus by the 
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unitary invariance of Choi-Kraus decompositions, it follows that EiPc — aiUPc for 
some scalars Ui. Hence we have A = (ctiaj) = \^p){'^p\, where is the vector state 
with coordinates a^, and so S'($,C) = S{A) = 0. 

On the other hand, suppose A — is rank one. Let y be a scalar unitary 

that diagonalises A. It follows that V induces a unitary change of representation 
for $ via ^ from {Ei} to {Fj}. But since FAT^^' is diagonal, only one Fj, say F, 
is non-zero, and hence by unitary invariance we have FPc = UPc for some unitary 
U. Thus $ o T^c = W o "Pc, and the result follows. ■ 

From an operational perspective, the numerical value of the entropy of a code 
allows us to quantify the number of ancilla qubits needed to perform a recovery 
operation. Specifically, the Choi rank, D, of a map $ gives the minimum num- 
ber of Kraus operators necessary to describe the map or, by Stinespring's dilation 
theorem |34j . the dimension of the ancilla required to implement $ as a unitary. 
The rank of A then gives the number of Kraus operators necessary to describe the 
action of the map restricted to C and thus the number of Kraus operators, M, 
necessary for a recovery operation in the usual measurement cum reversal picture 
of recovery. Again by Stinespring's dilation theorem we need an M-dimensional 
ancilla to implement the recovery as a unitary, and thus this requires logAf qubits. 

Hence the entropy is equal to zero for a unitarily correctable code, in which the 
action of the noise is unitary and thus requires no ancilla to implement the recovery 
operation. If the code entropy is positive, then any state of the code can potentially 
evolve to any one of multiple locations in the system Hilbert space under the action 
of the noise <i>. This fact has to be compensated by the recovery operation 4*. The 
maximal entropy case for a particular <i> is characterised by evolution to each of 
these locations (M = D by Lemma[T]) with equal probability (by Theorem[3]). Here 
the entropy and thus the number of qubits in the ancilla will be logD. 



2.4. Stabilizer Example. As an example from the stabilizer formalism, consider a 
three-qubit system with the usual notation Xi, Zi, i — 1,2,3, for Pauli operators [1]. 
The singie-qubit stabilizer code with generators {Z1Z2, Z2Z3}, is spanned by |0l) = 
|000) and ~ The set of operators {1, Xi, X2, X3} form a correctable 

set of errors for this stabilizer thus we can consider a channel, a three-qubit bit- 
flip channel, comprised of these operators - for example, the channel with Kraus 

p-q~r)l, E2 = ^Xi, E3 ^ 

'irXg. Using Pc = |000)(000| + |111>(111|, Q tells us that the error correction 
matrix is 

/ l{3-p-q-r) \ 


ig 

\ ir 



operators Ei 



hqX2 and E4 — 



A = 



The entropy of this code is therefore 



5(A) = -i(3-p-g-r)logi(3-p-g-r) 



-p log -P ~ -g log -g - -r log -r. 



As we would expect, the minimum entropy is achieved when p = q = r = 0. 
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The maximum entropy occurs when p = q = r = j, which we might also expect 
since this puts the auxihary density matrix A into the maximally mixed state. Here, 
the Choi rank of our map is 4 as is the rank of A, and the spectrum of A is equally 
balanced. So C is non-degenerate for $. Further, as we expect from Theorem[3l the 
entropy is log 4 = 2. Since the rank of A is the number of Kraus operators required 
for a recovery operation we find that in order to dilate the recovery operation to a 
unitary process, we need a 4-dimensional space, or equivalently a 2-qubit ancilla, 
which matches the value of the entropy. 

We next consider a variant of this noise model, in which the third qubit undergoes 
no error and the first and second qubits are flipped with probability equal to that 
of no error. Thus the map $ has noise operators {/, Xi,X2}, each weighted by 
The entropies for a trio of correctable codes for $ are given in Table [1] (vectors are 
assumed to be normalised). The first code is non-degenerate, and thus yields the 
maximal entropy for this noise model. The second code has positive entropy less 
than the maximum since it is partially degenerate for $. Indeed, one can check that 
Xi (and /) act trivially on the code, whereas X2 maps the code to an orthogonal 
subspace. The final code shows zero entropy since it is fully degenerate for $. In 
fact, as can be directly verified it is a decoherence-free subspace for 

Table 1. Codes for the noise model $ = -^{I,Xi,X2}- 



Qubit Code C = {|0l) 


,|1l)} 


Entropy S{^,C) 


{|000),|111)} 




log 3 


{|000) + |100),|011) + 


|111>} 


logs - 1 


{|000) + |100) -1- 1010) - 


Hiio), 





|011) + |111) + |001) + 


1101)} 





2.5. Entropy of a Subsystem Code. We next consider an extension of code 
entropy to the case of operator quantum error correcting subsystem codes. We 
shall only introduce this notion here and leave a potential further investigation 
for elsewhere. Subsystem codes were formally introduced under the umbrella of 
operator quantum error correction [31 [T^] , a framework that unifies the active and 
passive approaches to quantum error correction. These codes now play a central 
role in fault tolerant quantum computing. 

Let H he a Hilbert space of finite dimension N. Any decomposition H — {AiS) 
B)®JC determines subsystems A and B oiH. Given a quantum operation $ acting 
on H, we say that a subsystem B oi H is correctable for $ if there exists maps ^ 
on H and ta on A such that 

(12) *o $075^5 = (ta® ids) o'Pab, 

where Vab{p) = PabrPab and Pab is the projection onto the subspace A® B. 

Subsystem codes generalize standard (subspace) codes ^ in the sense that sub- 
spaces may be regarded as subsystems with trivial ancilla (dim A = 1). Thus, it is 
natural to suggest that a notion of entropy for subsystem codes should generalize 
the subspace definition. We find motivation for such a notion through the main 
result of [4] alluded to above. To every correctable subsystem for there are sub- 
systems C and B' ^ B, a map tc^a from A to C and a unitary map Vb'\b from B 
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to B' such that 

(13) <^oVab^{tc\a®Vb'\b)°Vab, 

Recall that the maximal output entropy of a channel ^' is given by S{^) := 
maxp S{'^{p)). This motivates the following. 

Definition 4. Let $ be a quantum operation with correctable subsystem B that 
satisfies (IT2t . Then we define the entropy of B relative to $ as the maximal output 
entropy of the associated ancilla channel, S{t(j^a) from p^ . 

Observe that this generalizes Definition [21 since a density operator may be re- 
garded as a channel from a one-dimensional input space. The minimal entropy 
case is characterized by the ancilla channel tc\a having range supported on a one- 
dimensional subspace. Such codes are unitarily correctable, in fact as subspaces, 
but the converse is not true. Instead, in the more general subsystem setting, the 
minimal entropy case is described by the associated ancilla subsystem A undergo- 
ing "cooling" to a fixed state. In principle one should be able to conduct a deeper 
analysis of subsystem code entropy. We leave this as an open investigation for 
elsewhere. 



3. Entropy of a Code for Binary Unitary Channels 
A binary unitary channel has the form 

(14) p' = $(p) (1 - p)WipWl + pW2pWl 

where Wi and W2 denote two arbitrary unitary operators and the probability p 
belongs to [0, 1]. 

It is clear that the problem of finding an error correcting code subspace C for 
the above map is equivalent to the case 

(15) p" ^<^u{p)^{l~p)p + pUpU^ 

where U = W\W2- The number M of Kraus operators is equal to 2, with Ei = 
\/l — pi and E2 = y/pU ■ Thus the error correction matrix A is of size two and 
reads 



(16) A = 



1-p -p)A 



where A is a solution of the compression problem for U 

(17) PcUPc = XPc ■ 

The set of solutions to this problem can be phrased in terms of the higher-rank 
numerical range of the matrix U . The rank-/c numerical range of U is defined as 

(18) nk{U) = {XeC\ PUP = XP for some rank-fc projection P} . 

Given a dimension-A: that defines the size of the desired correctable code, each A 
in flk{U) corresponds to a particular correctable code defined by the associated 
projection P that solves pT)) . The following is a straightforward application of ([7]). 

Proposition 5. Given a binary unitary channel $, there exists a rank-k correctable 
code for $ if and only if the rank-k numerical range of U is non-empty. 
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Figure 1. Higher rank numerical range flk for unitary matrices 
U describing a bi-unitary channel: two qubit system, a) example 
2 with A = 0, b) case with r = |A| >> for which code entropy is 
smaller, c) two qutrit case with A G ^^{U) chosen to maximize its 
modulus r and to minimize the code entropy S{A). 




Thus, the problem of finding the correctable codes for a given binary unitary 
channel can be reduced to the problem of finding the higher-rank numerical range 
of C7. This problem has recently been solved in its entirety [Te lfTTlfTSlfTQl 12011211122] . 
Most succinctly, in terms of the eigenvalues (j{U) of U, the kth numerical range of 
U is the convex subset of the unit disk given by 

(19) nkiU) = fl conv(r), 

rCcr{U); |r|=Ar-fc+l 

where conv {Ai, . . . , Am} is the set of linear combinations A = tiAi + . . . + tm^m 
such that J2Y=i = 1 ^^'^ ^ 0- Figures l.a and l.b depict the case of a generic 
two-qubit unitary {N — 4) with k = 2, while Figure l.c shows the case of a generic 
two-qutrit unitary noise {N = 3x3 = 9) with A: = 3. 

With a particular A in-hand, straightforward algebra provides us with the spec- 
trum of the matrix (jl6p . 

(20) A± = i(l±v/l-4p(l-p)(l-|AP)), 

which allows us to calculate the entropy of the code, 

(21) S'($,C) = 5(A) = -A+logA+ - A_logA_. 
We are then led to the following: 

Theorem 6. A minimum entropy rank-k code for a binary unitary channel $(p) = 
{l—p)p+pU pU^ corresponds to any code for which the magnitude |A| of the compres- 
sion values A € Qk(U) is closest to unity, while the maximum entropy corresponds 
to \X\ closest to zero. Moreover, a code with minimal entropy can be constructively 
obtained. 

Proof. The first statement follows directly from an application of first and second 
derivative tests on (|2ip . constrained to the unit disk. A minimal entropy code 
can be explicitly constructed based on the analysis of higher-rank numerical ranges 

[M [HI Ha Hi 1201 USUI]. ■ 

Example 7. As an illustration of the code construction in the simplest possible 
case {N = 4, A: = 2), let U be the unitary with spectrum depicted in Figure l.a, 
a{U) — {zk — exp(fc^) : k — 1,3,5,7}. Let j^'fe) be the associated eigenstates. 
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f/jV'fc) — ^fcl^Afc)- In this case we have r22(t/) — {0} so A = 0, and one can check 
directly that a single qubit correctable code for $ is given by C span{|(/ii), |02)}, 
where 

102) = ^(1^-2) + IV'4)) ■ 

For a concrete example, in the case that p = 0.01, ([2T|) yields a code entropy of 
S'($,C) = 0.081. 

The general case requires a more delicate construction, nevertheless it can be 
done. The "eigenstate grouping" procedure used above can be applied whenever fc 
divides iV. For instance, in the generic = 9 and fc = 3 case depicted in Figure l.c, 
a single qutrit code can be constructed for all A in the region 0.-^(1]). The states 
|</)i), i = 1,2,3, can be constructed in an analogous manner by grouping the nine 
eigenstates for \J into three groups of three, and writing A in three different ways 
as a linear combination of the associated unimodular eigenvalues Zj, j = 1, . . . , 9. 

However, without going into the details of this construction we can still analyze 
the corresponding code entropies. For simplicity assume the nine eigenvalues are 
distributed evenly around the unit circle with z\ = 0. By Theorem |6] we know 
that the entropy will be minimized for any A that gives the minimum distance from 
^z{U) to the unit circle. An elementary calculation shows that one such A, given 
by the intersection of the lines through the first and seventh, and sixth and ninth 
eigenvalues (counting counterclockwise) is approximately Aq — 0.092— 0.524 i. With 
the probability p = 0.01, the corresponding error correction matrix A has spectrum 
{0.007,0.993}. Thus, (UHl) and (HT)) yield the minimal qutrit code entropy for this 
channel as 

min S'($,C) = min 5(A) = ^(AA^An) = 0.060. 

dimC=3 Ae03((7) 

On the other hand, as A = belongs to Vl^{U), by Theorem [6] and (fT6|) we also see 
the maximal entropy for p — 0.01 occurs for any code with A = 0. In such cases we 
have A spectrum {0.01,0.99}, and hence the maximal entropy is 5(Aa=o) — 0.081. 

Changing focus briefly, if we fix an arbitrary unitary [/, then we could consider 
the family of channels determined by varying the probability p. It follows from ((2T|) 
that the channel with the correctable codes of maximal entropy corresponds to the 
p = \ channel, and the channels whose correctable codes possess minimal entropy 
correspond to p = and p = 1. Indeed, the value of A depends on U but not p, thus 
a given A will solve (jl7p for any p and so A can be chosen independently using the 
above theorem. The result once again follows from application of first and second 
derivative tests with p between and 1. 

The following results show that the entropy of a code for a binary unitary channel 
can be regarded as a measure of how close the code is to a decoherence-free subspace. 

Lemma 8. //$ is a binary unitary channel, then the sets of unitarily correctable 
subspaces and decoherence-free subspaces coincide. 

Proof. As proved in [J], for a bistochastic (unital) map a condition satisfied 
by every binary unitary channel, the unitarily correctable codes (respectively the 
decoherence-free subspaces) for $ are imbedded in the fixed point algebra for the 
map o $ (respectively $), where is the Hilbert-Schmidt dual map of 
In particular, it follows from this fact that the former set is given by the set of 
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operators that commute with U and W , whereas the latter is the set of operators 
that commute with U. By the Spectral Theorem these two sets coincide. ■ 

Theorem 9. Let ^ be a binary unitary channel. Then there is a rank-k code C of 
zero entropy, S'($,C) = 0, for^ if and only if there is a k- dimensional decoherence- 
free subspace for $ if and only if there exists X G flk{U) D (^{U). 

Proof. A fc-dimensional decoherence-free subspace for $ corresponds to an eigen- 
value A of [/ with multiplicity at least k (see [?] and the references therein); that 
is, A G rife(J7) n cr{U). The rest follows from the lemma and previous theorem. ■ 

In order to further illustrate these results, consider again the case of an arbitrary 
two-qubit system [N ~ 4). The correctable codes with largest entropy are those 
with p ~ \ and so the spectrum of A reads 

(22) A± = i(l±|A|). 

In the two-qubit case, the complex number A is given by the point inside the unit 
circle at which two diagonals of the quadrangle formed by the spectrum of U cross 
(see Figures l.a and l.b). Consider a special case of the problem where U has a 
doubly degenerated eigenvalue, so that |A| = 1. For example, U could be any (non- 
identity) element of the two-qubit Pauli group. Then the spectrum of A consists of 
{1, 0} which implies 5(A) — (despite p having been chosen for the largest entropy 
correctable codes). Hence A is pure and there exists a decoherence free subspace - 
the one spanned by the degenerated eigenvalues of U . 

In general, for binary unitary channels one may use the entropy (|2ip as a measure 
quantifying to what extend a given error correction code is close to a decoherence- 
free subspace. For instance, any channel (|15p acting on a two qubit system and 
described by unitary matrix U = W1W2 of size 4 may be characterized by the radius 
r = |A| of the point in which two diagonals of the quadrangle of the spectrum cross. 
The larger r, the smaller entropy S'($,C), and the closer the error correction code 
is to a decoherence-free space. 

The code entropy can be also used to classify codes designed for a binary unitary 
channel acting on larger systems. For instance in the case of two qutrits, iV = 3x3 = 
9, one can find a subspace supported on fc = 3 dimensional subspace. The solution 
is by far not unique and can be parametrized by complex numbers A belonging to 
an intersection of 3 triangles, which forms a convex set of a positive measure. From 
this set one can thus select a concrete solution providing a code C, such that r = |A| 
is the largest, which implies that the code entropy, S'(<I>,C), is the smallest - see 
Figure I.e. Such an error correction code is distinguished by being as close to a 
decoherence-free subspace as possible. 

4. Conclusions 

We have investigated a notion of entropy for quantum error correcting codes 
and quantum operations. The entropy has multiple natural realisations through 
fundamental results in the theory of quantum error correction. We showed how the 
extremal cases are characterised by unitarily correctable codes and decoherence- 
free subspaces on the one hand, and the non-degenerate case determined by the 
Choi matrix of the map on the other. We considered examples from the stabilizer 
formalism, and conducted a detailed analysis in the illustrative case of binary uni- 
tary channels. Recently developed techniques on higher-rank numerical ranges were 
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used to give a complete geometrical description of code entropies for binary unitary 
channels; in particular, the structure of these subsets of the complex plane can 
be used to visually determine how close a code is to a decoherence-free subspace. 
We also introduced an extension of code entropy to subsystem codes, and left a 
deeper investigation of this notion for elsewhere. It could be interesting to explore 
further applications of the code entropy in quantum error correction. For instance, 
although quantum error correction codes were originally designed for models of 
discrete time evolution in the form of a quantum operation, generalizations to the 
case of continuous evolution in time [3511361157] have been investigated. Further, we 
have investigated perfect correction codes only, for which the error recovery opera- 
tion brings the quantum state corrupted by the noise back to the initial state with 
fidelity equal one. Such perfect correction codes may be treated as a special case 
of more general approximate error correction codes [33l Ell El] . Another recent 
investigation [lU] includes analysis that suggests the measurement component of 
recovery may prove to be problematic in quantum error correction, and hence may 
motivate further investigation of unitarily correctable codes. 
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